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We consider the theory of feedback control of a Bose-Einstein condensate (BEC) confined in a 
harmonic trap under a continuous measurement constructed via non-destructive imaging. A filtering 
theory approach is used to derive a stochastic master equation (SME) for the system from a general 
Hamiltonian based upon system-bath coupling. Numerical solutions for this SME in the limit of a 
single atom show that the final steady state energy is dependent upon the measurement strength, 
the ratio of photon kinetic energy to atomic kinetic energy, and the feedback strength. Simulations 
indicate that for a weak measurement strength, feedback can be used to overcome heating introduced 
by the scattering of light, thereby allowing the atom to be driven towards the ground state. 
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I. INTRODUCTION 

The atom laser is the most coherent source for atom 
optical experiments [HE]. However, the applicability of 
the atom laser as a tool for fundamental research is lim- 
ited by noise that broadens the linewidth. This noise 
is due, in part, to instability in the spatial mode of the 
Bose-Einstein condensate (BEC) from which the atom 
laser beam is outcoupled. Excitations of the BEC spa- 
tial mode exist under general preparation conditions [3], 
and are generated when the atom laser is continuously 
pumped m |5] . One promising solution to this difficulty 
is to drive the BEC towards a stable spatial mode via the 
use of measurement feedback control. A control scheme 
based on a measurement of position has been shown to 
successfully cool a single atom in a harmonic trap close to 
the ground state from any initial state [6 . However, al- 
though this control scheme can be engineered by placing 
the atom in a cavity [7 , it is unclear whether it could be 
generalised to a many- atom BEC. In this paper we con- 
sider a feedback control scheme based upon dispersive 
imaging, a technique that has already been implemented 
in multiple BEC laboratories. 

Feedback on an atom laser was first applied in a single 
mode model to reduce phase noise caused by the interac- 
tions [H [9] . Improving the modal stability of a BEC using 
feedback control was first examined using a semiclassical 
model where it was shown that the system could be 
stabilised with feedback by modification of the trapping 
potential and the condensate's nonlinearity. This feed- 
back scheme was then applied to a more realistic model 
of an atom laser that included pumping, damping and 
outcoupling [5 . While the semiclassical approximation 
allowed for an examination of the multimode behaviour 
of the atom laser, the effect of coupling the system to a 
measuring device was ignored. The measurement backac- 



tion was included in a model of a trapped single particle 
considered by Doherty and Jacobs [7 , who showed how 
a position measurement arises from placing an atom in 
a cavity, and solved the optimal control problem for an 
initial Gaussian state. It was later shown that the filter 
equation could be solved, and the atom could be cooled 
from any arbitrary state |6j. The position measurement 
relied upon the assumption that the atom is trapped in a 
region small compared to the wavelength of light within 
the cavity. However this assumption is false for a mod- 
estly sized BEC in an optical cavity. Moreover, even if 
this type of position measurement could be engineered 
for a BEC in theory, either by situating the BEC in a 
cavity or otherwise, this does not imply that it would be 
easy to implement in practice. Indeed, a weak position 
measurement of a condensate has not been experimen- 
tally realised. There is thus a clear preference towards 
developing a feedback scheme that uses a well-established 
technique of measurement. 

There are two commonly used techniques for measur- 
ing a BEC of alkali atoms. The first method, termed 
absorption imaging ^ , shines near-resonant laser light 
on the condensate. Those photons which interact with 
atoms from the condensate are absorbed, leaving a 
'shadow' which can be detected using an array of CCD 
cameras. Thus absorption imaging gives a measurement 
of the column density (the number density integrated 
along the line-of-sight of the laser) of the BEC, which 
gives information about the spatial distribution of atoms 
in the condensate. The key advantage of absorption 
imaging is that any data obtained is independent of the 
intensity of the light, the time of exposure and many 
properties of the CCD array. A downside, however, is 
that the absorption of photons heats the atoms suffi- 
ciently to destroy the BEC. 

Phase- contrast imaging [TTJ [12] is an alternative 
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method of imaging that uses hght highly detuned from 
resonance. The interaction of the BEC with the hght 
gives a phase profile for the condensate, which can be 
used to reconstruct the density profile. In principle, 
the detuning can be sufficiently large such that a mea- 
surement of the density is minimally destructive. The 
theoretical limits for the signal generated from this dis- 
persive measurement scheme are not in fact different 
from absorption imaging in the case of low optical depth 
p!3| iMj [15]. However, for a BEC, it is far easier to 
achieve near optimal sensitivity for a given level of spon- 
taneous emission in a dispersive measurement. Exper- 
iments showing successive repeated measurements with 
phase-contrast imaging [11] suggest that it may be possi- 
ble to use this measurement technique non-destructively 
on the timescale needed to perform the feedback control. 
In this paper we show that measurements of a BEC via 
the technique of phase-contrast imaging can be used to 
construct a feedback control scheme. 

In Sec. |TI| we present our model of the system and 
the associated stochastic master equation (SME) for the 
quantum filter. After adiabatically eliminating the ex- 
cited state we recapture the SME considered by Dalvit et 
al. in the limit where the size of the condensate is much 
larger than the wavelength of light [16 . The form of 
feedback, which is modelled by the inclusion of a control 
Hamiltonian, is also outlined. In Sec. Ill the quantum 



filter is numerically solved in the limit where the atomic 
sample contains only a single atom. The dependence 
of the atom's final steady state energy on the measure- 
ment strength, ratio of photon recoil energy to typical 
atomic kinetic energy and form of feedback is discussed. 
Attempts at a numerical solution for the quantum filter 
under the semiclassical limit are also elucidated. 



II. MODEL AND FILTER DERIVATION 
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FIG. 1: (Colour online) Diagram illustrating how one could 
perform a non-destructive density measurement on a BEC. 
The BEC is illuminated with highly off-resonant laser light. 
The interaction between the field and the atoms is registered 
as a phase shift on any light scattered from the condensate. 
Such a phase shift is measurable by homodyne detection in 
the phase quadrature. 



The system under analysis is a BEC magnetically con- 
fined in a harmonic trap and illuminated with an off- 
resonant coherent field directed along the z-direction (see 
Fig. [T]). Light scattered from the condensate is detected 
by an array of homodyne detectors. The total Hamilto- 
nian for the combined system is 



sys 



B + Hi . 



(1) 



For simplicity, it is assumed that the atoms have been 
configured such that a transition between only two states 
is possible. These levels will be labelled the ground {g) 
and excited (e) states. Under this assumption, the sys- 
tem Hamiltonian can be written as 



i7,vs= / d^r^l{r)Ha{r)^g{r) 

+ / d^r^l{r){Ha{r)^hw^)^e{r), (2) 



where uoq is the resonant frequency of the atoms and '^i(r) 
is the field operator that annihilates a boson from the 
zth atomic level at position r, obeying the commutation 
relation ['^^(r), '0j(r')] = 5ij5^{r—r'). Ha{r) is the single 
particle Hamiltonian 
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For simplicity, we have taken the interaction energy be- 
tween atoms to be negligible, i^controi(^) is the single par- 
ticle control Hamiltonian, which we will specify shortly. 
The Hamiltonian for the electromagnetic field is 

Hb = j d^phc\p\ [dl{p)d,{p) + dl{p)d,^{p)) , (4) 

where e{p) and s'{p) are the two unit vectors required to 
describe the polarisation of the field, p = (kx^ky^kz)^ and 
[di{p)^ d^j{p')] = 5ij5^{p—p'). It can be assumed that the 
wavelength of the driving laser is much greater than the 
Bohr radius of an atom. Thus the interaction between the 
electric field and an atom can be adequately modelled by 
approximating the atom as a dipole. In this case the most 
general interaction Hamiltonian Hj for a field of two-level 
bosonic atoms interacting with the electromagnetic field 
is ill, 



Hi 



Jd'rxPlir) 
- J d'ri^lir) 



where dij is the transition electric dipole vector between 
the ith and jth level of the atom, defined as 



dij 



(6) 



where \i) denotes the state vector for the ith. level. We 
choose the phase of the internal states such that deg is 
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real, which imphes that d^g = dg^. E is the quantised 
electric field operator, which can be expressed in terms 
of a basis of creation and annihilation operators as [18] 



E{r) =^jd'pY.sl {-Me'"-"- - h-c] e{p). 

(7) 

Importantly, the coupling of the system to the bath 
allows us to make a number of standard quantum op- 
tical and reservoir approximations, thereby leading to 
tractable equations. We demonstrate in appendix [A] 
that making such approximations yields the following Ito 
quantum stochastic differential equation (QSDE) for the 
unitary of the system and bath: 



If we perform a homodyne measurement of the phase 
quadrature on the laser light after it has interacted with 
the atoms, then it can be shown that the best-estimate 
(in the least-squares sense) tt^ of any system observable 
X is given by the equation [20l [2T] 



d(j{t) 



J dx (^M{x)dA\x,t)- M\x)dA{x,t)^'^U{t), 

(8) 



dT:t{X) = TTt{L[X])dt -i j dx {i^t (xM{x) - M\x)X^ 
-TTt{X)nt (m{x) - M\x)) } dW{x, t), (11) 
where we have defined the Lindblad generator 
C[X] =^-[H,X] + j dx {m\x)XM{x) 

-\{M\x)M{x),X}\. (12) 



where 
H 



= J dx{i^lix)Haix)i^gix)+i^l{x) {Ha{x) + hA)Mx)} 
^[dge • €(0)] J dx {i^l{x)4^e{x) " h.C.) , 



dW{x, t) is the classical Wiener increment (i.e. Gaussian 
white noise). It satisfies dW{x^t)dW{x' ^t) = 5{x — x')dt. 
This noise is the random wavefunction collapse that cor- 
rupts the homodyne measurement signal. By defining the 
conditional density operator pc by 7rt(X) = Tr[Xpc] we 
can construct the following stochastic master equation 
(SME) for the quantum filter: 
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(9) 



and 



dpc = --[H,pc]dt^ J dxV[-iM{x)]pcdt 

+ / dxH[-iM{x)]pcdW{x,t), (13) 



M{x)= dx^l{x)^e{x)V{x-x), 
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for detuning A = loq — col and photon fiux of the laser 
Fq. zq is the characteristic length scale of the conden- 
sate in the tightly trapped z and y directions (i.e. zq = 
yjhjmujz)' The quantum Wiener increment dA{x^t) 
obeys the Ito differential rule dA{x, t)dA\x' ^ t) — 5'^{x — 
x')dt [H]. This term models the noise introduced into 
the system due to vacuum fiuctuations in the bath. The 
inner product between the dipole and polarization vec- 
tors can be shown to equal ^9] 



V[c]p = cpc^ - ^{c^c.p}, 
n[c]p = cp^ pc^ - Tr[(c + c^)p]p, 
for any arbitrary operator c. 

A. Adiabatic elimination 

For the trapped BEG under consideration the detuning 
of the laser is necessarily large. This allows us to adia- 
batically eliminate the excited state. We first transform 
to the following dimensionless units: 



Xo 



Ki — 



t' = tA, 



(14) 
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[p-dgef/\p\'- 



(10) 



The above QSDE corresponds to a 'cigar'-shaped con- 
densate where the BEG is tightly confined in the z and 
y dimensions. We present the one-dimensional version 
of the unitary evolution simply because one dimensional 
simulations require less comput national power. In ap- 
pendix [a] a more general two-dimensional derivation for 
a 'pancake' shaped condensate is given. 



where xq = ^JhJmuJr is the characteristic length scale of 
the trap in the x direction, which has trapping frequency 
cjt- The field operators are made dimensionless by the 
transformation V^i(x) / ^f^- Define the following 

parameters: 
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where dge = \dge\. Then Hamiltonian ([9l) and conditional 



master equation ( 13 ) can be written as 



H 



(15) 



For large detuning, any atom excited by the laser spends 
a relatively short amount of time in the excited state be- 
fore returning to the ground state. Moreover, there are 
very few atoms in the excited state in comparison to the 
ground state. Hence, on longer timescales it will appear 
that the population of excited atoms is tiny and changes 
very little. Thus, on this slower timescale we can approx- 
imate d^lje/df ^ 0. After making this approximation, 
eq. (18) gives 



and 



dpc = -i[H,pc]dT' + ^ / diV[M{0]pcdT' 



MO 



(19) 



y 

d^n[M{^)]pcdW{^, r ), (16) biguity upon substituting eq. ([19 



respectively, where 

HaiO = Ha{x)/nWT 

M{i) = -iJ de^PliOMOr.ii-C), 



f,(0 = 




7(^x) = / ^ dny^-r2 / 



for w = Note that the Wiener increment has 

been rescaled by mapping dW{x^t) dW{(,, r')/\^xoA. 
This preserves the ^-correlation dW{(,,r^)dW{(,\r') = 
S{£^ — £,')dr' . We interpret Fgp as the rate at which a 
single atom spontaneously emits into the bath, and Q 
as the Rabi frequency for the atomic system. We as- 
sume that the detuning A of the laser is much larger 
than other characteristic frequencies in the system, i.e. 
A Q^Tsp^ooT' Furthermore, the intensity of the laser 
is sufficiently large such that ft Fsp,cJT- 

We begin the adiabatic elimination by calculating the 
evolution of ^/^g, which can be found from eq. (11) by 
using the Heisenberg equation: 



dr' 



(17) 



Note that the terms proportional to Fgp /A and y^Fgp/A 
are small and have thus been neglected. Furthermore, 
those terms in H proportional to Ut/A can be ignored, 
as they are small compared with the terms proportional 
to Q/A and unity. Thus 



df 



,lot 



MO, / deMO—MC) 



.uot 



MO,i— J dc {i'i{e)Me) - h.c?j 

o . 

(18) 



As i/jI and iljg do not commute, there is an ordering am- 

into H. However, only 
one possible ordering yields a valid master equation: 

H = I d^i'liOHoiOMO + ^ I d^i'liOHoiOMO 
d^i^liOHaiOMO, (20) 



where the term proportional to 1/A^ is very small, and 
has hence been ignored. The conditional master equation 
simplifies to 



dpc = —i[H, Pc]dr + a 



I dmMaio] 



pcdr 



+ V^J d^n[MampcdW{^,T), (21) 

where we have chosen the more convenient time scahng 
r = (jjTt and 



M 
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± Sp ^ £, 
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(22) 



It is now clear that there are two key dimensionless 
parameters upon which the system depends. The pa- 
rameter a represents the strength of the measurement. 
For a larger a, more information is obtained in a fixed 
time. However, there is also more measurement backac- 
tion due to the linear scaling with a of the decoherence. 
T] is the Lamb-Dicke parameter, which is clear when it is 
written as = 27rxo/A = hko/po^ where A is the wave- 
length of light emitted from the atoms and po = Vhuom 
is the characteristic momentum spread in the harmonic 
trap. Thus 77 is proportional to the relative momentum 
'kick' an atom gets from scattering a photon. 

As a final note, typical experiments of the kind de- 
scribed above operate in the regime where zq ^ \ (i.e. 
w ^ 1). Thus the exponential exp[— tL'(l — nzY]i which 
appears in the integrand in the definition of 'y{nx)-, is 
tightly peaked about = 0. Physically, this is 

indicative of the fact that photons are strongly scattered 
in the z direction, and therefore have very little momen- 
tum in the x and y directions. Hence we are justified in 
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expanding = y 1 — (^^ + second order within 

the exponent: 

0.4 . (23) 



(^x + ^y) i^x + ^2/) 

* ' " —2k?: 



^0 



in which case (tZz — ko)'^ ^ (^^ + ^^)^/4^o- Furthermore, 

we can approximate tZz ^ 1 and [dge - sip)] ^ <^^e with 
very httle effect on the form of 7(/^a?). Finahy, we can 
extend the hmits of integration over hZy (and, incidentahy, 
those for integral over hz^ which defines r(^)) to ±oo, to 
give a much simpler functional form for 7: 



dtzy e-^N exp [-w{4 + i^lf/^] . (24) 



In this regime our result is in agreement with the master 
equation considered in [16]. 



B. Control 

Up to this point we have left the form of control com- 
pletely arbitrary. Indeed this is the advantage of sepa- 
rating the control problem into estimation and control 
stages. Since the filter provides a best estimate of the 
system state, these control terms can be a function of 
any system observable. Thus one is free to concentrate 
on the design of an effective feedback scheme, secure in 
the knowledge that the filter for the system is indepen- 
dent of the choice of feedback. We consider a general 
control Hamiltonian of the form 



J^co„trol(0 = E""W^"- 



(25) 



n=l 



Haine et al. [4] performed a semiclassical analysis that 
showed that the change in energy of the BEC was always 
non-negative for the choice 



(26) 



dt ""2 

for positive constants c^. In this paper we will be pri- 
marily concerned with a feedback control consisting of 
only the first term in eq. (25). More precisely 



HcontroliO = {P) C (27) 

A feedback control of this form was used in [6^ and [7 . 
Physically, this control represents an adjustment of the 
trap minimum such that the motion of the atom is damp- 
ened. Such a control could be implemented experimen- 
tally via the use of changing magnetic fields. 



III. SIMULATION 

A. Single atom limit 

Important insight into the behaviour of a BEC under 
the above mentioned control scheme can be gathered by 



considering the physical limit of a single atom. In this 
case the density operator can be written as 

where the coefficients ^(^, are given by 
The evolution of these coefficients is given by 

de{i,a = {modpi'Ho\o)- (28) 



Substituting eq. (13) into eq. ( |28[ ) yields the following 
SME: 



dQc 



dr -\- a J dhix 



V 



Qcdr 
(29) 



where 



Ha 



cZH^(/t^,T) 



3a 

2^ 
1 



Ci (p) X 



die-''^-^dW{i,T). 



dW^Kx-,^) is the Fourier transform of the Wiener incre- 
ment. It behaves somewhat differently to the traditional 
Wiener process: 

dW"" {Kx-,^)dW{K'^,T) = 5{kx — i^'x)dT 
dW{nx,T)dW{i<i'^,T) = 5{nx + i^'x)dT. 

It is possible to express dW{nx^ r) in terms of the Wiener 
increment dW^nx^r). In particular, 

dW{^x.T) = ^ - 1) {dW{^x.T) + idW{-Kx, r)) . (30) 

Ultimately, we are interested in how the average en- 
ergy of the system E \{x^) + (p^)] /2 varies over time; 
in particular the steady state value for the energy. In 
the following analysis we judge the effectiveness of the 
control based upon how close the system is cooled to 
the ground state energy TujOt/^ and the time taken for 
the system to reach a steady-state. The energy was cal- 
culated by finding a numerical solution to the stochas- 
tic Schrodinger equation corresponding to eq. (29). The 
numerical integration was performed by using the open 
source software package xpdeint, which is a new version 
of the xmds package [22 . 

Simulations revealed three important features of the 
system. The first relates to the measurement strength 
parameter a. For a sufficiently large a the final average 
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FIG. 2: (Color online) Plots of the energy for simulations of 
500 paths. Simulations were for a normalized Gaussian func- 
tion centered at x — 2.0. Parameters chosen were ci = 2.0, 
w = 3000.0, and 77 = 6.0. The measurement strength (in order 
from bottom to top) is a = 2.0 (red) , 10.0 (blue), 20.0 (green), 
40.0 (black), 80.0 (cobalt), 160.0 (magenta). Fuh hues repre- 
sent the mean, while dotted lines indicate the standard error. 
Notice that those simulations with a < 20.0 converge to the 
same average steady-state energy. 



FIG. 3: (Color online) 500 path simulations of the energy for 
an initial normalized Gaussian state centered at x = 2.0, it; = 
3000.0, a — 10.0, and ci = 2.0. The Lamb-Dicke parameter 
(in order from bottom to top) is 77 = 4.0 (red), 6.0 (blue), 
8.0 (green), 10.0 (black). Full lines represent the mean, while 
dotted lines indicate the standard error. This plot shows that 
relatively small increases in 77 result in large increases in the 
final average energy. 



steady state energy increases with increasing measure- 
ment strength. This is because the measurement has a 
greater backaction on the atom for a larger measurement 
strength, which corresponds to an increased heating rate. 
Thus, a small a is required for a low steady-state en- 
ergy. A caveat, however, to choosing a small a is that 
less information is obtained about the system per unit 
time. This translates to an increase in the time it takes 
for the energy to reach a steady-state. For optimal con- 
trol, one would like to balance these considerations by 
choosing an a that cools close to the ground state energy 
on a timescale much smaller than that of the experiment. 
However, an additional constraint is that the cooling does 
not continue to get better as a decreases. There is a 
threshold value ac, where for any a < ac the final aver- 
age energy is the same as that for ac- Indeed, the only 
effect of decreasing a lower than ac is that it takes longer 
before the energy reaches a steady-state. A demonstra- 
tion of this dependence on a is shown in Figure [2] If no 
other constraints are taken into consideration, then this 
suggests that ac is in fact the optimal value for a. 

The second point to consider is the dependence on the 
final steady state energy on the Lamb-Dicke parameter 
77. For small 77, the energy imparted to the system dur- 
ing the scattering of the light is negligible. However, 
for large rj the centre of mass dynamics of the atom is 
greatly influenced. This introduces additional heating 
into the system. It was found that for larger values of 77, 



the higher the final average energy (see Figure [s]). Fur- 
thermore, while the final average energy scales roughly 
linearly with a, it seems that it scales at a greater rate for 
77. Preliminary investigations seem to indicate that this 
scaling is quadratic. However, a more thorough study is 
required in order to draw a more accurate conclusion on 
this scaling. For typical trap frequencies (cc; ~ 0.1 Hz - 1 
kHz), we can see that for a rubidium atom, rj will take 
values somewhere between 1 and 400. These results, cou- 
pled with the limited ability to reduce the height of the 
plateau by decreasing a, indicate that this measurement 
and control scheme may only effectively cool a trapped 
atom for strong trapping potentials. 

The third interesting feature involves a consideration 
of the effect of adding additional control terms of the 
form ( [25| ) to the Hamiltonian. As demonstrated in Fig- 
ure [4j additional feedback terms increase the effectiveness 
of the control. In fact this is only true for larger values 
of 7^; if 7^ < 1 then the higher order controls have little 
effect on the final average energy. This is precisely be- 
cause for a large 77 there is significant coupling between 
the different modes of the atom. This coupling is due to 
nonlinearities introduced by the trapped atom experienc- 
ing different electric fields at different positions in space. 
Hence energy can be removed from those higher order 
modes which are not directly affected by the feedback 
control. This is in contrast to a situation with small 
77, where higher order modes remain unaffected by the 
control of the lower order modes. Thus, despite the ad- 
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FIG. 4: (Color online) 100 path simulations of the en- 
ergy in the single particle approximation with feedback terms 
proportional to x (red/highest), x and x^ (blue/middle), 
and (green/lowest). The parameters chosen were 

w = 3000.0, a = 2.0, and 77 = 8.0. Full lines represent the 
mean, while dotted lines indicate the standard error. The 
initial states chosen were normalized Gaussian functions cen- 
tered at X = 2.0. The constants of proportionality for each 
feedback term were set to 2.0. 



ditional heating due to a large 77, it may be possible to 
effectively cool the atom close to the ground state with 
the introduction of additional control terms. 



B. Mean-field approximation 

Although the single particle limit has its uses for qual- 
itative understanding, multi-particle effects are impor- 
tant in many BEG phenomena. Indeed, the dynamics 
resulting from collisions between the atoms simply can- 
not be included in a single particle model. The next step 
towards modelling the above situation for a BEG with 
interactions is to consider the quantum filter (13) under 
the mean-field approximation 



This yields the following Ito equation of motion for the 
mean- field wavefunction (/)(^): 



■ a- 



(31) 



In contrast to the single atom limit, numerical simula- 
tions of eq. (31) are not convergent. We conjecture that 



the difficulty stems from an inherent mismatch between 
the mean-field approximation and the type of measure- 
ment associated with phase-contrast imaging. This imag- 
ing gives information about the average density of the 
condensate, and is in effect a measurement of the num- 
ber of atoms in the BEG. Thus over time, the continuous 
measurement projects the BEG state towards a number 
state. However, upon making the mean-field approxima- 
tion it is assumed that the number variance is always that 
of a coherent state. This incongruity with the measure- 
ment may explain why simulations of the quantum filter 
under this semiclassical approximation give non-physical 
results. 



IV. SUMMARY AND CONCLUSIONS 

) 

This paper has investigated the prospect of generating 
a stable spatial mode for a trapped BEG by using feed- 
back control to cool the condensate close to the ground 
state. In particular, the state of the BEG was estimated 
using a realistic measurement scheme based upon non- 
destructive imaging. We presented the full Hamiltonian 
for the system and bath, and derived a stochastic master 
equation for the conditional density matrix of the system. 

This SME was simulated in the single atom limit, and 
it was shown that a feedback scheme which adjusted 
the trap minimum would bring the atom to a steady- 
state energy. However, the precise value for this final 
energy and the time taken to bring the atom to that 
energy depends upon the measurement strength, Lamb- 
Dicke parameter (ratio of scattered photon's kinetic en- 
ergy to atom's kinetic energy) and feedback strength. 
Simulations revealed that decreasing the measurement 
strength decreased the average steady state energy, but 
only up to a point. Below a certain critical value the 
steady state did not decrease; only the time taken to 
reach a steady state increased. A more problematic is- 
sue is the additional heating introduced into the system 
for values of the Lamb-Dicke parameter corresponding to 
some weakly trapped BEG experiments. However, it was 
demonstrated that much of this heating could be offset by 
the introduction of additional control terms proportional 
to higher powers of x. 

The derived SME was also simulated under the mean- 
field approximation, for a condensate with no interaction 
energy. These simulations did not converge. We hypothe- 
sise that this lack of convergence is due the measurement 
projecting the state to an eigenstate that is no longer 
coherent, which is a dynamical effect at odds with the 
approximation itself. We are currently investigating the 
behaviour of this filter under a number-conserving semi- 
classical approximation. Ultimately, however, a definite 
answer to this question may only be obtained via a full 
field calculation. Such full field calculations may be pos- 
sible using stochastic techniques that are also currently 
under investigation [23 . 
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APPENDIX A: DERIVATION OF ITO 
STOCHASTIC DIFFERENTIAL EQUATION FOR 
UNITARY EVOLUTION 

In this appendix we present a derivation of eq. ([8| from 
the total Hamiltonian ([T]). The key physical approxima- 
tion involved is that the optical field acts as a Markovian 
reservoir. This requires a modal restriction of the atomic 
field, which is equivalent to allowing the possibility of a 
light pulse to travel away from the BEC and have no fur- 
ther interaction with it. Although it is possible to make 
reservoir approximations for the light when a finite num- 
ber of modes all interact with the same reservoir, we shall 
simplify our geometry by assuming that the magnetic po- 
tential is sufficiently tight in the z direction such that the 
BEC is highly restricted in this dimension. This gives a 
'pancake'-shaped condensate, and restricts the BEC to 
the occupation of a single mode z direction. That is 

i,g{T)^g{z)i,g{x), 4(r-)«5(^)e''=°^4N, (Al) 

where x = (x^y). '^i{x) is interpreted as the two dimen- 
sional field operator for the zth state at position x. The 
inclusion of the factor e*^°^ represents the phase shift an 
atom receives upon being excited by a photon of momen- 
tum on the order of /cq = ooq/ c. Under this approximation 



X (^2Pl{x)2Pg{x)e-'^^' ^h.c^ . 



(A2) 



It is reasonable to assume that the density profile of the 
condensate in the z direction is a normalised Gaussian of 
width zq. Specifically 



\9iz)\' 



1 



(A3) 



Notice that before the dimensionality of the conden- 
sate was restricted there was one mode of light per mode 
of the condensate. After restriction, there are two in- 
finite dimensions of radiation for every mode of BEC. 
This allows us to treat the z component of the electro- 
magnetic field as a reservoir. This is most simply done 
by first transforming the integral over into an integral 
over frequency uj. The frequency is related to the wave 
number in the z direction by the relationship 



(A4) 



for fe = {kj;,ky). This change of variables gives 

[dge ■ £{k,ui)] 



Hi 



■oo / 



: (^L\k)d{k,uj)G{k^ - ko) ^ L{k)d{k,u)G{k^ + k^) 
-L\k)d\k,u)G{k^ + ko) - L{k)d\k,u)G{k^ - ko)^ 



where 



and 



m = I 



ih-x 



(A5) 



(A6) 



G{k,) 



1 



/27r 
exp 



dz\g{z)\ e 



2 ^ — ikoz ^ — ikzZ 



(A7) 



The annihilation operators have been undergone the 
rescaling d{p) y^^l^^Jco d{k,Lj) to ensure that the 
commutation relation 

[d{k,u;),d'^{k\u;')] =S\k- k')8{u - u') (A8) 

is preserved. Note that even though ^ikz give the same 
cj, the operators d{k^kz) and d{k^—kz) act on different 
Hilbert spaces. However, given most of the scattering 
will be in the positive z direction, we have ignored the 
contribution due to d{k,—kz). The Hamiltonian for the 
electromagnetic field becomes 

Hb = j d^k J duhujd\k,uj)d(k,uj). (A9) 

We move into the interaction picture with the unitary 
transformation 



Ui{t) = exp \^Hb + hcjL J d'^xi)l{x)ij^e{x)'j , 

(AlO) 

where is the optical frequency of the laser. This trans- 
forms the operators as follows: 

a±{k,oj) a±(fc,w)e-''^*, L{x) L(a;)e-''"^*. 

(All) 

Now the coupling of the atoms to the electromagnetic 
field will occur predominantly in a narrow frequency 
range cjq — < co < cjq -\- for some <^ coq. The 
rotating wave approximation can thus be used to neglect 
those terms which are rotating quickly. Hence 



Hi 



d^k 



dujK,{k^uj) (L^ {x)d{k^L 
n JujQ-e ^ 



OL)t 



-L{x)d/ {k^ uo)e^ 



(A12) 
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where we have defined the strength of couphng between 
the system and photon bath as 



dge ■ e{k,Lo) G{kz - ko)- 



2(27r)3c2fc^fe:o 

(A13) 

The domain of integration over k has been restricted to 
Q = {k : c\k\ < ujo} to ensure that uj is never complex. It 
is assumed that the coupling strength is roughly constant 
in frequency space around the resonant frequency. That 
is, K.{k,uj) ^ hz{k,uJo). Furthermore, let 



— / 



dix) d(k^ uj)e' 



■ (A14) 



Then 

Hi = ih j (fkV2^K{k,LOo) (L(fc)a(^)1'(fc,t) - /i.e.) . 

(A15) 

So far no direct assumptions have been made about 
the nature of the electromagnetic field coupled to the sys- 
tem. This field can be well approximated as a classical 
light field with quantum vacuum fiuctuations. Make the 
replacement d^^\k^t) d^^\k^t) -\- f{k^t)^ and approx- 
imate the entire bath as a vacuum state. The interaction 
Hamiltonian becomes 

Hi = ihJ d^kV2^i^{k,uo) \^(^L{k)d^^^'^{k,t) - h.c.^ 

+ (L{k)r{k,t) - L\k)f{k,t))]. (A16) 

In the interaction picture, the evolution of the unitary 
operator [// is 



d 
It 



Ui{t) 



Hi 



(fx^l{x) {Ha{x) + hA)Mx) Ui{t), 



(A17) 



where A = uq — lul. Substituting (A16) into this expres- 
sion gives 



-U^^Ht) = {-lH^ [ d^kV2^K{k,u;o) 

X (^L{k)d^^^'^{k,t)-h.c}j}u^^\t), (A18) 

where 

H = I d^x^l{x)Ha{x)^g{x) 

+ J d^X^lix) {Ha{x)^hA)i^e{x) 

ih J d^kV2^n{k,ij) (^L{k)r{k,t) -h.c.y 

(A19) 



The superscript {0) has been used simply to highlight 
the dependence of the unitary on ^. It is reasonable to 
assume that the BEC is much smaller in the x and y 
directions than the spatial size of the coherent beam. 
This allows the laser to be adequately approximated as a 
plane wave propagating in the z direction with frequency 
col- Thus f{k^t) ^ \^ToS'^{k)^ where Fq is the photon 
number fiux of the laser. Under this approximation the 
Hamiltonian H reduces to eq. 
Now, note that as ^ ^ oo. 



d^^\k,t) d{k,t) 



1 



In vacuum (a(fc, t)) = and (^d^{k^ t)d{k, s)) = 5{t — 5), 
and can therefore be identified as quantum white noise. 
Given that the coupling is weak, in the sense that 
nik^ix}^) <C 1, eq. (A18) indicates that the timescale on 
which the system evolves will be slow. Since the fre- 
quency range of coupling is narrow, it can be assumed 
that ^ cjo)]^- Hence the system is well described 
by taking as practically infinite [18 . The importance 
of this observation is that the quantum analogue of the 
Wong-Sakai theorem states that in the limit of {k^t) 
approaching quantum noise (i.e. 9 00) the solution 
to eq. (A18) approaches the solution to the Ito quantum 



stochastic differential equation (QSDE) 



-H 



i j d^k-i(k)L\k)L(k)^ dt 



+ j d^k^^i(1^(^L{k)dA\k,t)-L\k)dA 

(A20) 



where we have defined 



7(Aj) = 27r / d^k[n{k,ujQ)f, (A21) 

and dA{k^ t) is the quantum Wiener increment. It satis- 
fies the property dA{k^ t)dA\k' ^ t) = S{k — k')dt. For an 
heuristic development of this theorem, see [20^. A more 
rigorous treatment can be found in [24l [25]. Thus, to 
a good approximation, the u nitarv for the total system 
is given by t he Q SDE (A20). Using the definition of L 
given in eq. (A6) and 



dA(x 



d^kA(k,t)e 



-ikx 



(A22) 



QSDE (A20) can be written in the form 

dU{t) = I ^-^iJ - i J (fx M\x)M{x)J dt 

+ J fx (^M{x)dA\x,t) - M\x)dA{x,t)'^^U{t), 

(A23) 
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where 



H 



CpX 'tpl{x)Ha{x)xpg{x) 



In this form, it is possible to integrate out the exphcit 
dependence on ky. Define 



7(fc^)=27r / dky[Ki{k,uJo)] , 



(A31) 



2(27r)3£o ' 



'-^[dge ■ e(0)] / d^x (i^l{x)Mx) - h.. 



Also consider the one dimensional quantum Wiener 
increment A{kx^t)^ which satisfies the relation 



and 



M{X)= j d^x'^{x'){ije{x')V{x-x'), 

V{x) = [ 



27r 



(A24) dA{kx,t)dA^k^^,t) = 5{kx - k'^)dt. So 

(^j d^kV2^ni{k,uoQ)dA{k,t)^ (^j d^k' V2^Ki{k\i^Q)dA{k\t)^ 
d^k2iT[i^i{k,ujQ)fdt 
dkj:-f{kx)dt 



J-} 



ko 



ko 



dkx y^j{kx)dA{kx,t) 



ko 



ko 



dk',^^)dA{k',,t) 



for dA{x,t)dA\x.,t) = S^{x - x')dt. 

As a final note, a 'cigar'-shaped condensate can be con- This shows that eq. ( |A2Q[ ) obeys the same statistics as, 
sidered by assuming the BEC is tightly confined in both and is therefore equivalent to, 
the z and y directions. In this case 

^g{r) ^ g{z)h{y)^g{x), ^ g{z)h{y)e'^'^ipe{x), 



(A25) 

where tpi{x) is the one dimensional field operator for the 
ith state at position x and 



dU{t) = 

rko 

ko 



r^o 



2 



If 

- J dkx-f{kx)P{kx)L{kx) ] dt 



dkr 



\Hy)\' = 



-yVz'o 



(A26) 



For simplicity we have assumed that the density profiles 
in the y and z directions are both normalised Gaussians 
of identical width zq. A similar argument to that given 
above leads to the evolution of the one dimensional uni- 
tary 



^/^(^L{k,)dAHk,,t)-L\k,)dA{k,,t))^ 

(A32) 

As before, this can be rewritten to give eq. ([8|. 



dU{t) 



^i^- ^^^i2/c2^[A^i(/c,c^o)]'^^fe)^fe)) dt 

d^kV2^ni(k,uoQ) (^L(k^)dA\k,t) 

-L\k,)dA{k,t))]u{t). (A27) 
where the coupling constant is now given as 



tvi{k^ uj) 
for 

and 



2{2'Kfc^k,heo 



dge-s{k,oj) G{k^-ko)H{ky), 
(A28) 



H{ky) 



dy\g{y)\ e 



L{k,) = J dxi^l{x)Mx)e-"'''''. 



(A29) 



(A30) 
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